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ABSTRACT 



In the case of adiabatic flow of an ideal electron or ion gas in the 
presence of gravitational and electromagnetic fields, it is possible to 
derive a system of scalar equations of motion for the gas that provides 
an alternative to the Euler formulation. The new dynamical variables 
that replace the velocity are three scalars, one of which is a generaliza- 
tion of the flow potential and satisfies a generalized Hamilton-Jacobi 
equation. The other two are constants of motion that specify the intrinsic 
vc ticity, which is that part of the total vorticity that is not produced by 
the external fields, but rather is simply a residue of the initial conditions 
of the gas. It is this intrinsic vorticity that satisfies the generalization 
of the simple vorticity conservation theorem. The same formalism can 
be adapted to the case of constant-temperature flow, and to the flow of an 
arbitrary barotropic fluid. 



SCALAR FORMULATION OF IDEAL CHARGED GAS FLOW 


I. INTRODUCTION 

Euler's equation for an ideal electrically uncharged fluid, which is a vector equa- 
tion, can be reduced to a scalar equation in the case of potential flow. The conditions 
that must be fulfilled in order for potential flow to exist are the following: (1) The 
pressure must be a function of the fluid density alone. (2) The external force acting 
on the fluid must be the gradient of a scalar potential. (3) The curl of the fluid velocity 
must everywhere vanish. 

According to the Helmholtz (or Kelvin) vorticity conservation theorem, if the third 
condition is satisfied at any given instant of time, it will continue to hold for all time. 

When the above conditions are satisfied, the velocity can be expressed as the 
gradient of a flow potential : 




( 1 . 1 ) 


The most important special case for which the first of the above conditions is 
satisfied is that of isentropic flow , for which the specific entropy is everywhere con- 
stant for all time. For this case, Euler's equation becomes 

- * S/*t = G + f.-v-jv 1 } (i.2) 


1 



where 0 is the potential for the external force (which throughout will be taken to be 
gravitation), and ft, is the specific enthalpy. 

In the case of time -independent flow, -*v&t reduces to a constant E having the 
dimensions of energy/unit mass, and (1.2) becomes 

E = G + *& + *$r v* ^ (i.3) 

which is Bernoulli's equation. 


It is interesting to note that 

DS — (i.4) 

where 

D— + X- £ < l - 5 > 

is the substantial time derivative. But from single-particle Hamilton- Jacobi theory we 
have the following relations: 

V.i = (1.6a, 

and 


Dl= = L ? 


(1.6b) 
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* where ^ is Hamilton's Principal Function, m is the particle mass, and L is the particle 
Lagrangian. Comparison of equations (1.6) with (1.1) and (1.4) indicates that the flow 
potential S plays a role analogous toj/^. Comparison of the potential energy terms 
in (1.4) and (1.6b) indicates that, as far as dynamical effects are concerned, the enthalpy 
plays the role of a thermal potential which is to be treated in the same way as the 
gravitational potential <5 . 

This analogy between the formulation of potential flow in terms of the flow potential 

S and the single -particle Hamilton- Jacobi formalism is further reinforced by using 

2 

(1.1) to eliminate V in (1.2). The resulting nonlinear particle differential equation for 
S has the form of the Hamilton- Jacobi equation: 

- * s /&t = G + < + ^ (vS)*. a-’) 

This equation provides a scalar alternative to Euler's equation in the case of ir- 
rotational, isentropic flow. Once (1.7) has been solved (for given (5 and 'ft), V,can be 
found from (1.1) and the fluid particle density ^ can be found from the continuity equa- 
tion as follows: 

D^e = - - V*S . (i.8) 

It would be desirable to generalize this scalar formulation in three respects: (1) by 
replacing the isentropic requirement that the entropy be everywhere constant by the 
weaker adiabatic requirement that the entropy be constant along any given flow line, 
in which case the pressure will be a function of the entropy as well as the density; (2) 
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by including the non-potential electromagnetic forces that act on an electrically charged 
fluid; (3) by including rotational as well as irrotational flow. This generalization will 
be carried out in the work that follows. We shall see, moreover, that a slight modifica- 
tion of the formalism suffices to include the constant temperature case for which the 
temperature must be constant on a given flow line (but may vary from line to line), and 
the case of an arbitrary barotropic fluid for which the pressure is an arbitrary function 
of the pressure alone. 

Throughout the derivation we shall consider only an ideal charged fluid. In the 
case of an electron gas, for example, this means that we shall neglect not only its self- 
viscosity, but also the frictional drag between it and the ion gas which necessarily oc- 
cupies the same region of space. 

H. EULER FORMULATION 

Our starting point is the familiar Euler equation for an ideal charged fluid, such 
as an electron or ion gas in which viscosity and heat conduction are ignored, in the 
presence of a gravitational potential (r and electric and magnetic fields £ and B : 

m D = - m VG- - f ' 7 P + + 0»/c) ^x < 21 ) 

where m and ^ are the particle mass and charge, P the pressure, and f the particle 
density. 

This equation must be supplemented by the continuity equation: 
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or 


+ v-(ev)= o 


Ce.za) 


\)JL* f s— VV . (2.2b) 

AM 

These two equations must be further supplemented by an equation that specifies 
the nature of the thermodynamic constraint on the flow. In the case of isentropic flow 
of a perfect gas for which the ratio of specific heats is V, this thermodynamic side- 
condition is 



Constant , (Isentropic Flow) 


( 2 . 3 ) 


For the weaker requirement that the flow be adiabatic, i.e. that the specific entropy be 
constant along a flow-line, the thermodynamic equation is 

D^p = -yoA,e =-*v.v . (24 

(Adiabatic Flow) 


A necessary first step in deriving a scalar alternative to Euler's equation is to 
eliminate f from (2.1). This is accomplished by means of the following thermodynamic 
identity: 

f 1 V P = m T V> — m ( 2 . 5 ) 

where T is the temperature, and A* and are the specific entropy and enthalpy 

respectively. 
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We shall also eliminate the electric and magnetic field intensities from (2.1) in 

o 

favor of the vector and scalar potentials A and A by means of the familiar relations 


E = - v a° - *A/ C *t 


(2.6a) 


and 




(2.6b) 


Using (2.5) and (2.6), (2.1) can be cast into the following form: 

It me A) “ " + * + ^A°) 

)+3rc&] . 


(2.7) 


The Thermal Function 

It is possible to absorb the termTVAin (2.7) into the other three terms by elimi- 
nating the temperature T in favor of a scalar function 3*", called the thermal function, 
that is defined by the following equation: 


DJ= T . 


( 2 . 8 ) 


The utility of the thermal function follows from the identity 

(o^kva) = (DA)(vr) - 

(2.9) 

+ V(^^/5t) + vx(2^xVj') , 

We now formally introduce the requirement that the flow be adiabatic: 

= O , (2.10) 
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* Using (2.8) and (2.10) in (2.9), we arrive at the desired expression forTVA.: 

TV> = 

( 2 . 11 ) 

+ i l '(P X ^) . 

Using this expression to eliminate the termTVA' in (2.7), we arrive at the following 
form of Euler’s equation: 

= -vc<T+*+iv%iA*-A*J7at) 

( 2 . 12 ) 

+ ^x[(vxv>4-i e S.+(^xv3 r )] _ 

Before pursuing the formal development further, it would be well to discuss the 
physical significance of the thermal function^*". From (2.8) we note that ^could be re- 
garded as the reading of a clock that moves with the fluid and whose rate is proportional 
to the temperature X. 

A further insight into the physical meaning of 2 an be gained from a comparison 
of (2.8) with (1.6b). The latter equation states that Hamilton's Principal Function £ 
can be regarded as the reading of a clock whose rate is the Lagrangian L which, in the 
absence of any potential energy, is just equal to the kinetic energy. But the temperature 
T is just equal to (except for a proportionality constant) the kinetic energy of random 
thermal motion. This suggests that plays a role analogous toX» except that it refers 
to the energy of random thermal motion rather than to the macroscopic kinetic energy. 
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Referring to (1.6a), we see that this analogy suggests that should in some way 
be related to a thermal momentum flux in the fluid, i.e. the momentum associated with 
a heat flux. This thought is reinforced by writing (2.12) in the following form: 

*£At = +- V,X (VX£) <213<0 

where 

P - + (Vc)A 4* (2.i3b) 

and 

•f = *r\ + -v ^ v i + (%/U)A° — _ (2.13c) 

We note that has the form of a generalized canonical particle momentum in which the 
term mA7^is the thermal contribution. This fact, together with the above-mentioned 
analogy of withX , suggests that the thermal function could perhaps be used to 
describe heat flux as well as temperature (i.e. heat energy density). It is evident that 
(2.8) does not completely determine 'S since it does not specify the difference in the 
values of S associated with different flow lines at a given instant of time. Because 
we have postulated an ideal fluid, and so are ignoring heat flux, we shall simply take 
the point of view that the degrees of freedom in that are not specified by (2.8) may be 
arbitrarily specified according to mathematical convenience, i.e. they may be regarded 
as inconsequential constants of integration. 
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Thermal Gauge Invariance 


Part of the indeterminacy in S' can be described in terms of a gauge indeterminacy 
that is completely analogous to what we encounter in the case of the electromagnetic 
scalar and vector potentials. We note from (2.13b) and (2.13c) that if we add an arbitrary 
function of to S', then |> is increased by the gradient of some scalar function ty', and 
jP is diminished by . Thus 


and 



+• v if 


(2.14a) 


(2.14b) 


But if and IP satisfy Eiiler's equation (2.13a), then so do and r. Thus Eider's 

/VW 

equation does not suffice to specify the gauge function V • As in the case of the electro- 
magnetic gauge, the choice may be made on the basis of mathematical convenience. We 
shall later see that a particular choice suggests itself in a very natural way. 


Intrinsic Vorticity 

Inasmuch as the definition of vorticity for an uncharged isentropic fluid is = 

(VX^vy,^ , a natural generalization would consist of replacing the particle momen- 
tum mv with the canonical momentum £ defined in (2.13b). Doing this, (2.13a) becomes 

^t*At = - wyf + x* <2-15*) 
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where 


u> =. (v x Vzm . (2.15b) 

^ <%*> 

We can gain an insight into the physical significance of w by noting that (2.15b) 
can be written in the following form: 


A =n t + 


(2.16a) 


where 




(2.16b) 



(2.16c) 


A* = ± (vTx v» . (2.i6d) 

A is the local angular velocity of the fluid. JQl l , whose magnitude is the Larmor fre- 
quency, represents the contribution to A that results from the diamagnetic response of 
the charged fluid to the imposition of the magnetic field JJ. A^is an analogous contri- 
bution to A that is thermal, rather than magnetic, in origin. It vanishes in the isentropic 
case - q), and even for the weaker adiabatic requirement, the undetermined 
degrees of freedom of 3 * can often be used to make A^vanish. Thus to is that part of 
fl that is not actively produced either by the magnetic field or by thermal effects, but 
rather is the part of A that is intrinsic to the fluid and is to be associated with the 
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initial conditions of the fluid. For this reason it will be called the intrinsic vorticity . 
(The definition (2.15b) indicates that the name canonical vorticity would also be 
appropriate.) In order to distinguish between to andl^, the latter will be referred to as 
the rotation, and the term vorticity will be reserved for tu. 

The importance of u> arises from the fact that it, rather than/L, obeys the gener- 
alized vorticity conservation theorem, which we shall now derive. 


Generalized Vorticity Conservation 
Taking the curl of (2.15a), and using (2.15b), we find 

- V X (*XU) . (2.17) 


It follows directly from (2.15b) that 


V * 0-> = O (2.18) 

/VS /WV # 


It is well-known that a vector which satisfies equations of the form (2.17) and (2.18) 
must also satisfy the following relation: 


°I 


u/.dS - O 


(2.19) 


where the integration is taken over an arbitrary surface S that is frozen into the fluid 
and moves with it. 


11 



Equation (2.19) is the desired generalization of the vorticity conservation theorem. 
We note that for an uncharged isentropic fluid B - Va = O, so co reduces to Jl and (2.19) 
becomes the familiar Helmholtz (or Kelvin) theorem. 

m. VORTICITY CONSTANTS OF MOTION 

As a preliminary to deriving a more convenient formulation of the vorticity con- 
servation theorem, we express the generalized canonical momentum in terms of 

( 2 ) 

three scalar functions $ ,M, and<p : 

P = V ^ n (3.1) 

An ^ ' 

where for later convenience S has been multiplied by the particle mass *n. Referring 
to (2.14a) we see that wS ihcfv<»es the indeterminate thermal gauge function. The sig- 
nificance of M and (f> becomes apparent when we substitute (3.1) into (2.15b): 

to = (tm)' 1 (V M) . (3.2) 

Thus M and $ determine the intrinsic vorticity QJ . We must now find the equations of 
motion satisfied byM and . These two scalar equations will then replace the vector 
equation (2.17). (Note that, because a> must be solenoidal, it can have only two degrees 
of freedom.) 

To arrive at the desired equations of motion for M and <J) , we substitute (3.1) into 
(2.15a): 
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(2>MAt)v<J> - 0><t>/dt) VfA - - V (f>+ m ^S/^t + M^to/bt) 

(3.3) 

+2m vxw . 

A%A AH. 

We now make a specific choice of thermal gauge as determined by the equation 

= f/m + (M/tv») 

(3.4) 

— (%/m) A —XbX/'b t + (M/m) ^ <t>/ • 


For this choice of thermal gauge (3.3) becomes 

V(t> - (b<t>/bt) VM - Zm v^X uj . (3.5) 

But from (3.2) we have 


= (l-V»)VM -(v-VM)V() . ( 3 . 6 ) 

Subtracting (3.6) from (3.5) we find 

(dm) V<|> - (D<t)VM = o . (3.7) 

Taking the cross-product of this equation with VM, and using (3.2) we have 

(PM) (*> = O . (3.8) 


Thus, if Co o, 


DN\ - O 


(3.9) 
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Similarly, 


D <t> = o 


(3.10) 


If co = o , then |> must be the gradient of a scalar, and IA and <b may be set equal to 
arbitrary constants (or zero) with no loss of generality. 

Equations (3.9) and (3.10) are the desired scalar equations of motion that describe 
the vorticity and replace the vector equation (2.17). These equations show that, for the 
choice of thermal gauge made in (3.4), KA and cj> are constants of motion of the fluid. 

It is this fact that makes this choice of gauge the most natural one. 

From (3.2) it is evident that the product M<#> must have the dimensions of angular 
momentum. With no loss of generality we may assert that M has the dimensions of 
angular momentum and that <P is dimensionless. Thus M may be regarded as an angular 
momentum, in some way related to the vorticity, that is conserved along a flow line. 

We may think of <f) as the initial value at t = O of one of the particle coordinates (in 
dimensionless form) , whose memory is retained by the particle for all time. 

As an illustration, consider the simple case of rigid cylindrical rotation of an 
uncharged fluid. The constant angular velocity is A ~fi k where k is the unit vector 
along the cylindrical axis. Let r be the distance from the axis and «p the azimuthal 
angle. Then if 

X. 

|V\ =: vv, ilr (3.lla) 
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and 

<$> - Cp-Jlt (3.11b) 

it is easily verified that (3.9) and (3.10) are satisfied. It follows from (3.2) and (3.11) 
that 

u> = tzm)' 1 (VM)X(V<t>) . (3.12) 

In this case the intrinsic vorticity <o and the rotation A are identical, and M is just the 
particle angular momentum. For a more general vortex, however, in which the mag- 
netic force were taken into account, we see from (2.16a) that For such a vortex, 

AH 

M would not be the total particle momentum, but only that part of it caused by the 
intrinsic vorticity, i.e. IA~>~yu~>r . 

Since in the present example the azimuthal angle Cp of a particle is given by 
cp=jt.t+<p 0 , where ip 0 is the particle's azimuthal coordinate at t = o , it is evident that 
<P = <p e . Thus, in this case it is the azimuthal coordinate whose initial value the par- 
ticle remembers. 


Finally, we note that for this example 


W7A — m/Lr U? x 

*** AH 


(3.13) 


where <p, is the unit vector in the azimuthal direction. Referring to (3.1), we see that 
MV<|> is the contribution to p that is to be associated with the intrinsic vorticity to , 

IS 

which in this case just equal to the angular velocity A • Since the right side of (3.13) 

A AV\ 
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is equal to the particle momentum which in this case is identical to , we see that 

"" <v»A 

for this example there is no potential contribution V(*v>S)to f> . 

Transition to w = O 

From (3.2) we note that if either M or (J> becomes everywhere constant, then u» 
vanishes. Moreover, because M and <P are constants along each flow line, if either 
one becomes everywhere constant at any instant of time (i.e. assumes the same value 
on all flow lines) , it must remain constant for all time. This corresponds to the fact 
that if, at any instant of time, oj * O in (2.17), then its time derivative also vanishes 
everywhere, and so uj must remain zero for all later time. 

The use of the scalar constants of motion M and $> to describe the vorticity pro- 
vides an insight into the manner in which the transition from ^ / o to«=0 takes 
place. In a fluid free of small scale instabilities such a transition could never come 
about. If, however, the Reynolds number is sufficiently large, turbulance mixes the 
fluid so that the values of the constants of motion M and <P (and/6/) on neighboring flow 
lines tend to equalize. This process could be pictured as a diffusion of the constants 
of motion from regions in which they have large values into regions in which the values 
are smaller. If the turbulent mixing and the time lapse since t-o are large enough, 
the fluid will have forgotten the various values of and <p that were assigned to each 
of the flow lines at t = O . In this case, it will have made the transition to values that 
are constants throughout the fluid, which means that the fluid has permanently settled 
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into the condition characterized by w-o. When the value of a, becomes everywhere 
constant, the fluid has made the transition from adiabatic to isentropic flow. 

Even when abundant turbulant mixing occurs, however, the transition to constant 
values of M , 4> , and & will not take place if an external constraint or boundary condition 
prevents it. For example, a rotating paddle wheel could inject vorticity (i.e. large 
values of N\ and $) into the fluid in such a way as to compensate for the decrease in M 
and 0 caused by diffusion (i.e. turbulent mixing). 

The constraint could also be geometric in nature as, for example, in the case of a 
fluid blowing radially outward from a central point (e.g. the solar wind). In this case 
as time elapses the ever increasing distance between what had been neighboring blobs 
of gas at t-o makes mixing ever more difficult. 

The above two constraints could occur in an ideal (i.e. frictionless) fluid. If we 
admit the possibility of friction, we have an additional mechanism that can prevent the 
transition to to -o- From (2.16a) we see that, to the extent that A* is negligible, 
to - o means that JT = J^/xmc . Thus, in the case of a fully ionized plasma consisting 
of electron and ion gases, the two gases must rotate in opposite directions if each is 
individually to fulfill the condition to =o . But the friction between the two counter- 
rotating gases, which is proportional to the electrical resistivity of the plasma, opposes 
this relative rotation and tends to make it vanish. Thus in this case friction works 
against the tendency of turbulent mixing to reduce u> to zero. 
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These examples illustrate that in many situations the transition to the condition 
*g~°, A> - constant does not occur because of constraints of one kind or another. If, 
however, no such constraints are present, and if one is not interested in a transient 
solution that is dependent on a particular set of initial conditions, then the particular 
solution characterized by =c has greater physical significance than any 

other solution that satisfies the same boundary conditions, in the sense that it repre- 
sents the "steady-state" condition to which the fluid will be brought by turbulent mixing. 
Such "steady-state" solutions should be of particular importance in astrophysical 
problems in which the Reynolds number, and hence the turbulent mixing, is generally 
extremely large. 


Potential Flow 

Referring to (3.1), we note that if either M or <|> is a constant (or if M is an 
arbitrary function of <t>), the term MV 0 can be expressed as the gradient of a scalar, 
and so can be absorbed into the term V£ws) . In such a case is just the gradient of a 
scalar and = o . Such flow will be designated adiabatic potential flow . Referring to 
(2.13b), we see that for such flow 

VS = v + (V«c) A +A- v S ' . (i'M) 

(Adiabatic Potential Flow) 

If now A' should become constant, then the tervaj^vS in (3.14) becomes the gradient 
of the scalar which can be absorbed into the potential function S . Such flow will 
be designated isentropic potential flow. It is characterized by the condition 
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(3.15) 


VS = v+(V~c)A . 

(Isentropic Potential Flow) 

In the above discussion we noted that, in the absence of any constraints that would 
prevent its realization, the flow characterized by the conditions to =0 and a = constant, 
which is just isentropic potential flow, represents the "steady-state" flow brought about 
by turbulent mixing. Taking the curl of (3.15), we can characterize such flow by the 
relation 

z = - (Vlwc) JB (3.16) 

(Isentropic Potential Flow) 

or^l= JX L » i* e * direction of the local angular velocity vector is opposed to the mag- 

netic field vector and has the magnitude of the Larmor frequency. 

It is of interest to note parenthetically that this condition could not be fulfilled by 
a constant-pressure charged fluid, because in this case there is no pressure gradient 
to help provide the necessary centripetal force, which must then be provided solely by 
the magnetic force. This requires that i.e. the angular velocity must equal 

the cyclotron frequency (ZA») rather than the Larmor frequency. Referring to (2.16a) 
we see that the condition Si- iil. (in the case of isentropic flow for which Jl . = o) re- 
quires that u» = u> u . Thus constant-pressure isentropic flow of a charged fluid is 
characterized by the presence of an intrinsic vorticity whose magnitude is the Larmor 
frequency. (Note that if aj = -JL U , then XL= o . Thus if rotation is present *v=JJl t > and 
if the rotation is zero uj = -JL u . In either case \<£ I - | .) 
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IV. SCALAR FORMULATION 
Adiabatic Flow 

Using (2.13b) to eliminate in (3.1), we obtain the following expression for 

-Vws 

V : 

1= SS - CVmc)A -AVT + (*/«) V 4> (4.1) 

^ ^ AM AM ’ y 

. z. 

which is the generalization of (1.1). Using (4.1) to eliminate ^ Y in (3.4), we arrive 
at the generalization of (1.7): 

— — 0*/^) — G •V’fc + (tyrw)A° 

(4.2) 

+ *k fe $ - cv~c^-^r -4* Wjfl 1 

which will be referred to as the generalized Hamilton -Jacobi (H-J) equation. 

From (3.4) and (4.1) we find , making use of (3.9) and (3.10), that 

[)S=i.V t ’-[(j+ , £ + C Vvnc)(cA° - V'jftJ)] (4.3) 

which is a very natural generalization of (1.4) inasmuch as, except for the factor 

— 'A. 

0 - cvef] ' which is unity in the non-relativistic limit in which we are working, 
c/\°— V' A) i- s J u st the contraction of the 4-vector- potential and the fluid 4-velocity, 

AM ' 

which means that the extra term in (4.3) that is missing from (1.4) is just the contri- 
bution of the electrostatic potential in the fluid rest -frame. 
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The density ^ can be found from the continuity equation. Using (4.1), we have 


D.&\ ^ = — * V • v 

AW 

- -V l s +(VW»<0£A + 

which is the generalization of (1.8). 

The independent variables of the problem are S ,M , <|> , A», and-oS". The generalized 
H-J equation (4.2) plus the side conditions P f/ ' = = o provide four of the five 

necessary equations. (Note that the v, that appears in the D operator is to be either 
eliminated by means of (4.1), or else treated within the framework of an iteration 
procedure as a given vector function.) 


To complete the system of equations, we need an equation for^Tl In order to derive 
such an equation, we must know the thermodynamic properties of the fluid. For the 
sake of illustration, let us consider the fluid to be a perfect gas for which the ratio of 
specific heats is — c p/c v • Because for a perfect gas f^=CpT, we have using 


( 2 . 8 ) 



(4.5) 


(Even in the case of an arbitrary fluid, ^ can always be expressed as a function of 
T= 0-5^ and a^, and so is not an independent variable.) For a perfect gas 


a> = ( k /«>) £<*-!)■%£ -JUt] * t * constant (4.6) 
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where k is the Boltzmann constant. From (4.6), the continuity equation (2.2b), and 
the adiabatic condition (2.10), we find 

D£ = - fear- 0 V-v]< } ( 4 . 7 > 

which can be written in terms of 3 ? by means of (4.5): 

dM” = - [cy-Ov- v] t)J" . <4.8) 

This is the desired equation for^T* which completes our system of five equations for the 
independent variables S ,N> ,<p , and^. 

In the chart given in Figure 1 a comparison is made between the Euler and the 
scalar formulations. In the Euler formulation we may regard the three components 
of ^ as the dynamical variables, and P and f* as the thermodynamic variables. The 
dynamical equations are the three components of the Euler equation; the equation for 
P is the adiabatic condition (2.4); and the equation for f* is the continuity equation 
(2.2b). 

In the scalar formulation^ is replaced by the three scalars S , and <f>, and the 
Euler equation is replaced by the H - J equation and the two side-conditions DM= *©. 
In the case of potential flow, a significant simplification occurs in the scalar formula- 
tion in that Wand <p drop out. No corresponding simplification occurs in the Euler 
formulation. 
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The thermodynamic variables in the scalar formulation are & and X which must 
satisfy the equations ° and (4.8), which are roughly comparable in difficulty with 

the equations for P and p in the Euler formulation. The density P does not enter into 
the scalar formulation at all. After the problem has been solved, and a> and X (and 
hence "R. ) are known functions, P can be found from (4.6) which, solved for f , becomes 

K — n it e <4.9) 

where N is a normalization constant. 

When the flow is isentropic rather than adiabatic, the only surviving thermodynamic 
variable is & Thus the single equation (4.7) replaces the two equations for A/ and X* 
that were needed for adiabatic flow. The only corresponding simplification that occurs 
in the Euler formulation is that the algebraic equation (2.3) replaces the differential 
equation (2.4). 

In the case of isentropic potential flow, the only surviving variables in the scalar 
formulation are S and ’ft, and the H- J equation assumes the following simplified form: 

-^S/bt = G +-K. + (V m ) A° + -f[v S - (^e)A~j (4.10) 

This has exactly the form of the single-particle H - J equation except for the presence 
of 

It is obvious that the advantage of the scalar formulation is greatest in the case of 
isentropic potential flow. In the more general case of adiabatic non-potential flow that 
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differs only slightly from isentropic potential flow, a simple perturbation solution is 
possible, using a solution of (4.10) as the zero-order solution. This procedure is 
facilitated by noting that, using 0<t>-O,DJ=T , and (4.1), (4.2) can be written 


— = G + 0& — + (%/rn)A? + 


where in this form of the generalized J equation it is the specific Gibbs function 
, rather than'R , that plays the roles of the thermal potential. The terms 
~ T * and-^Vxf- in (4.11) constitute the perturbation. Because both of 

these terms will be small for flow that is only slightly different from isentropic potential 
flow, the change in 5 will be small, and we may linearize the term 
in calculating this change in S . 


Constant -Temperature Flow 

The term constant-temperature flow will be used to designate flow that satisfies 
the condition 

DT = O (4.12) 

which replaces the adiabatic condition (2.10). Flow satisfying the stronger condition, 

T = constant, will be designated isothermal flow . 

The derivation of the scalar formulation in the case of constant -temperature flow 
is completely analogous to that for adiabatic flow. We can eliminate the density in 
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Euler's equation (2.1) by means of the thermodynamic identity 


-I 


? V P — m V a. -f- m^VT 

^ Q M 


(4.13) 


where 


^ (4.14) 

is the specific Gibbs function. In order to transform the term tnAflT thus introduced 
into Euler's equation, we define a function X by means of the equation 

D-i * ^ . (4.15) 

By means of an identity analogous to (2.9), Euler's equation can be written in the form 


=-»vf>+ v x<vx£) 


(4.16a) 


where now, in the constant -temperature case, 

p - mv + (% lc) A — wT vj, (4.i6b) 

and 

IP — +■ y 4* -Jr v *f (^/m)A° 4* T (4.16c) 

Equations (3.1) and (3.2) are valid for constant-temperature flow as well as for adiabatic - 
flow. Thus we see that any adiabatic flow equation can be converted into the corresponding 
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constant -temperature flow equation by replacing(« ; ^ T, and 3 ) by (jr -T, ^ 
and respectively and making no changes in S ; ^ or oj. 

The generalized H" J equation for constant -temperature flow is 

- - T ^ At - (**/«) = <3r+fr + <V"0 A * 

(4.17) 

+i[v5 -(^/me)A +TVJ, +<M/ m )V«Q^ 

This equation must be augmented by the three side conditions DT= P = O m To 

complete the system of equations we need an equation for & analogous to (4.8). This 
equation, like (4.8) j will depend on the thermodynamic properties of the fluid. As before, 
we shall consider a perfect gas. Since in this case C. p T » we have 

< 3' = -R-T*> ^T(Cp-O^) (4.18) 

and so 

P^=-TD l i (4.19) 

where we have used (4.12). It is possible to arrive at an independent expression for 
from the thermodynamic identity 

P DP ” m D<J- -V- m,4» PT (4.20) 

where in the second step we have again used (4.12). Since for a perfect gas P— k^T , 
(4.20) becomes 
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0 f = (kT/rn) OA»f = - (kT/m)V.V 


(4.21) 


where in the second step we have used the continuity equation (2.2b). Finally, from 
(4.19) and (4.21) we have 

D^*(k/»)V-V (4-22) 

which is the desired equation forX that is analogous to (4.8). Thus we have completed 
our system of equations for the three dynamical variables 5 ? and <J) and the two 
thermodynamic variables ~T and^ . The density ^ is not one of the simultaneous 
variables, but once the problem has been solved, ^ can be found from 


?=NT V(»-I) e -^ ( 4 . 23 ) 

which is just (4.9) with the factor Cp absorbed into the normalization constant 

N . The velocity V can be found by substituting (3.1) (which remains valid for constant- 
/w \ 

temperature flow) into (4.16b): 


v = vs-Wmc)A + Tvi + . ( 4 . 24 ) 

The form of the constant -temperature H - J equation that is analogous to (4.11) is 

-'bs/bt =&■*■■€.+ Cym)A° + ^[vS-(V m c)A ] 2 

^ -i Z (4.25) 
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where now plays the role of the thermal potential. This fact is especially 

favorable in the case of a perfect gas since then *£. = c p"*~ and it is not necessary to 
use the more complicated relation (4.18) to express the thermal potential (which in the 
case of (4.17) is <^) in terms of the thermodynamic variablesT and £ . For a perturba- 
tion calculation for nearly isothermal potential flow, however, we would write ^ +TA/ 
in (4.25) instead of and would regard!^— kfTvj.4 as the perturbation. 


For isothermal flow (T = constant) the terms involving^ in (4.17) can be absorbed 
into the terms involving 5 simply by replacing 5 with the new flow potential 5=5 . 


The isothermal H - J equation can also be derived directly from the adiabatic 
equation (4.2) simply by specifying that where T" c is the constant temperature. 

Then 'hS'/'bt = and V-3T " -O ( and so (4.2) becomes 

— — (x4- ^ 4- (%/m) A° 

(4.26) 

+ tk s - C Vrn«) A + <w/*0 v*]* 

(Similarly, the isentropic H- J equation can be derived from (4.17) by specifying 

where ^ is the constant specific entropy.) Using D<t> = 0 and the expression 
(4.24) for V (without the termT V& which has been absorbed into the term V 5 ), we 
arrive at the following alternative form of (4.26): 

2. 2L 

=■ Cr4- f + CVm)AV(r[S;S- (MnAAj - [CM/m) V <^| . (4.27) 
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In the case of isothermal flow, the only thermodynamic variable is <j. , so we would 
use (4.21) as our thermodynamic equation instead of (4.12) and (4.22). The most con- 
venient expression for £ is found by substituting A. ~ Cp - <^/t (from (4.18)) into 


(4.23) yielding 


_ m 




(4.28) 


Barotropic Fluid 


In the case of a barotropic fluid there exists a thermal potential function 1^-7^ (f ) 
that is a function of the density alone such that 


f'vP=v» 

M* AM \ • 


(4.29) 


This relation can be used to eliminate ^ from Euler’s equation (2.1). The derivation 
of the scalar formulation is identical to that for adiabatic flow except that "ft must be 
replaced by Ij and the terms involving ~T ; & and do not appear. Thus the scalar 
formalism is identical to that for isentropic flow except for the replacement of "f\ by?f . 
The independent variables are ^ and ^ . The three dynamical equations are 

the H- J equation and the conditions DM= D$> -O . The thermodynamic equation for 

is derived by solving the relation V[- for f in terms of ^ , and then substituting 

this into the continuity equation P jUy, £ = ~ V/X to arrive at the desired equation 
forP^ . 
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V. DISCUSSION 


The scalar formulation has an obvious advantage over the Euler formulation in 
the case of potential flow, since two of the three dynamical variables (M and ) drop 
out, whereas in the Euler formulation it is necessary to retain three variables, namely 
the three components of . In simple one-dimensional problems this advantage vanishes 
since v has only one degree of freedom, but for three-dimensional problems it can be 
very significant. 

Even in the case of non-potential flow, the fact that two of the three dynamical 
variables in the scalar formulation, M and <$> , are constants of motion should make the 
fluid behavior easier to visualize and perhaps easier to calculate. 

The scalar formulation must always be applied within the framework of a two-fluid 
picture of the plasma, i.e. the dynamical equations of the ion and electron gases must 
be worked out individually. (A third, uncharged fluid can be added if necessary.) This 
procedure is more exact than the usual perfectly-conducting-single-fluid model in that 
the inertial effects of the electric current (essentially the electron flux) are taken into 
account. Because of the assumption that the electron and ion gases are ideal, however, 
we are neglecting the dissipative part of the plasma resistivity. 

In the Euler formulation the intrinsic vorticity to is never brought into the open 
since it is regarded merely as a certain differential function of the components of v, , 
which are the basic variables. In keeping with this approach, all functions of the fluid 
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flow, such as current density and frictional forces, are expressed as functions of v^. 
In the scalar formulation w (via the scalars M and <f> ) plays a central role, and £ is 
reduced to the status of a certain differential function of S f M ^ and ♦ . In keeping 
with this point of view, it would be natural to express such things as current density 
and frictional forces as functions of and<£ . Such a procedure might suggest 

some physically meaningful approximations that would simplify the overall problem. 
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Footnotes 


1. See, for example, L. D. Landau and E. M. Lifshitz, ’’Fluid Mechanics" (Addison- 
Wesley Publ. Co., Reading, Mass., 1959), p. 19. 

2. Such a representation is always possible for a continuously differentiable vector 
function. See, for example, L. Brand, ’’Vector and Tensor Analysis" (Wiley & Sons, 

New York, 1947), p. 230. 

The introduction of the scalars S, M, and 0 in equation (3.1) represents 
a generalization of a procedure first used by Clebsch (Crelle’s Journal, (1857) 
and J56 (1359) ) for an uncharged barotropic fluid. This procedure is 
summarized by H. Lamb, ’’Hydrodynamics" (Cambridge University Press, Cambridge, 
1932) , 6th ed. , p. 248, and by H. Bateman, "Partial Differential Equations 
of Mathematical Physics" (Cambridge University Press, Cambridge, 1932), p. 164. 
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Figure Caption 


Fig. 1. Comparison of Euler and Scalar Formulations. 
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